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Solar system constraints on f(T) gravity 
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ABSTRACT 

We use recent observations from solar system orbital motions in order to constrain f(T) grav- 
ity. In particular, imposing a quadratic f(T) correction to the linear-in-T form, which is a 
good approximation for every realistic case, we extract the spherical solutions of the theory. 
Using these spherical solutions to describe the Sun's gravitational field, we use recently deter- 
mined supplementary advances of planetary perihelia, to infer upper bounds on the allowed 
f{T) corrections. We find that the maximal allowed divergence of the gravitational potential 
in f(T) gravity from that in the teleparallel equivalent of General Relativity is of the or- 
der of 6.2 x 10~ 10 , in the applicability region of our analysis. This is much smaller than the 
corresponding (significantly small too) divergence that is predicted from cosmological obser- 
vations, as expected. Such a tiny allowed divergence from the linear form should be taken into 
account in f(T) model building. 
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1 INTRODUCTION 

According to accumulating observations of different kinds, the uni- 
verse is now known to be accelerating (Riess et al. 1 1199 81) . Perl- 
mutter et al. dl999l) ). Such a feature led physicists to follow two 
directions in order to explain it. The first one is to introduce the 
concept of dark energy (see Copeland et al. ( feOOq) and Bamba 
et al. a2012al) and references therein) in the right-hand-side of the 
field equations of General Relativity, which could either be the sim- 
ple cosmological constant or various new exotic ingredients (Ra- 
tra & Peeb les dl988t) . Wette rich dl988h . Caldw ell J2002h . Nojiri & 
Odintsov ( hOoJT Feng et al. J2005h . Guo J2005h . Cai et al. fcoich ). 
The second direction is to modify the left-hand-side of the general 
relativistic field equations, that is to modify the gravitational the- 
ory itself, with the basic extended gravitational theories known as 
/(.R)-gravity (see (De Felice (2010)) and references therein). Such 
extended theories can present very interesting behaviors; the cor- 
responding cosmologies have been investigated in detail (Nojiri & 
Odintsov ( 120071) ). 

A new class of modified-gravity theories, which has recently 
received attention by the scientific communi ty, is t he f(T) gravity 
(Ferraro & Fio rini d2007|). Fer raro & Fiorini ( 1200 81) . Bengochea & 






Ferraro ( 120091) . Linder J20ld) ). It is an extension of the old idea 
of the "teleparallel" equivalent of General Relativity (TEGR), in 
which one uses the curvature-less Weitzenbock connection instead 
of the torsion-less Levi-Civita one, and where the dynami cal ob - 
jects are the four l inearly independent vierbeins (Einstein ( 1928), 
Unzicker & Case J2005I) . Hayashi & Shirafuji dl979l) . In TEGR, 
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the torsion tensor is formed solely from products of first derivatives 
of the tetrad. Then, the Lagrangian density T can be constructed 
from such a torsion tensor under the assumptions of invariance un- 
der general coordinate transformations, global Lorentz transforma- 
tions, and the parity operation, along with requiring the Lagrangian 
density to be second order in the torsion tensor (Hayashi & Shira- 
fuji dl979l) ). Thus, in f(T) gravity one generalizes the above TEGR 
formalism, making the Lagrangian density a function of T, similar 
to the f(R) extension of Einstein-Hilbert action. 

The above /(r)-gravitational framework proves to lead to in- 
teresting cosmological behavi or and it has gained much attention 
in the literature (Mvrzakulo v (1201 ID . Yerzhanov et al J2010I) . Wu 
& Yu bOlOal) . Chen e t al. <|201ll) . Wu & Yu .201 lah . Bamba et 
al. feOlOh . Dent e t al. J2011l) .~Zneng & Huan g l2oilh . Ba mba e t 
al. J201ll) . Yang bOllh . Zhang et al. ( 1201 ll) . Cai et al. ( 1201 ll) . 
Chattop adhyay & Debnath feOllh . Sharif & R anifeOllh . Wei et al. 
d2011al), Fer raro & Fiorini (201 la), Wei et al. ( 201 lb), Capozziello 
et al. fcOllI) . Wu & Yu fcoilbl) . Ba mba & Geng .201l[) Gen g 
et al. d201lh . Wei d2012t) . Geng et al .2012I) . Wu & Geng d201ll) . 
Bohme r et al J2012I) . Karami & Abdo lmaleki feOllh, At azadeh & 
Darabi J2012J) . F arajolla hi et al. d2012l) . Yang et al. J2012I). Karami 
& Abdol maleki J2012a|), Karam i & Abdolmaleki J2012bj) . Xu et 
al J2012I) . Bamba et a l. J2012bl) . Setare & Hound jo 120121) . Liu et 
al. d2012t) . Wu & Yu |2010bl) . Bengochea J201lh, Gonzalez et al. 
12), Daouda et al. feOllh. D aouda et al. 1 2012a]) Boehm er et al 
2011 ). Daouda et al . .2012bl) . Ferraro & Fio rini J2011ch . Wang 
( 20111). Miao et al. ( l201ll) . Wei et al J2012I) . Ferraro & Fiorini 
( 2011a) ). One question that straightforwardly arises is what are 
the qualitative and quantitative /(r)-modifications that are allowed. 
Although theoretically one has an enhanced freedom, in practice, in 
order to confront observations one expects a small divergence from 
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TEGR, that is from the linear-in-r c ase wh ich coincides wit h Gen- 
eral Relativity. Indeed, in Wu & Yu ( feOlObl) and Bengochea J201lh 
the authors used data from cosmological observations in order to 
constrain the model parameters of two well-studied /(r)-ansatzes, 
namely the power-law and the exponential one, and they found less 
than 1 % divergence from TEGR. 

In the present work we use latest observations of solar system 
orbital motions in order to constrain the parameters of f(T) gravity. 
In particular, after extracting the spherical solutions of the theory 
we use them to describe the Sun's gravitational field and then we 
can use data from planetary motion in order to infer an upper bound 
on the allowed divergence from the TEGR. The solar system analy- 
sis imposes significantly tighter bounds than the cosmological ones 
(Wu & Yu d2010bl) . Bengochea J201lh ) as expected, and in this 
completely different observational region it verifies too that the di- 
vergences from TEGR (and thus from General Relativity) are very 
small. 

The plan of the work is as follows: In section [2] we briefly 
review f(T) gravity and in section [3] we extract the spherical so- 
lutions assuming small corrections to the linear-in-r scenario. In 
section|4]we use data form solar system orbital motions and we im- 
pose constraints on the model parameters. Finally, in section [5] we 
summarize the obtained results. 



2 F(T) GRAVITY 

In this section we briefly review f(T) gravity. Our notation is as fol- 
lows: Greek indices fj, v,... and capital Latin indices A, B,... run over 
all coordinate and tangent space-time 0, 1,2, 3, while lower case 
Latin indices (from the middle of the alphabet) i,j,... and lower 
case Latin indices (from the beginning of the alphabet) a,b,... run 
over spatial and tangent space coordinates 1, 2, 3, respectively. 

As stated in the Introduction, the dynamical variable of 
"teleparallel" gravity, as well as of its f(T) extension, is the vier- 
bein field e^(.i^). This forms an orthonormal basis for the tangent 
space at each point xf 1 of the manifold, that is e^ • e B = r] AB , where 
J]ab = diag(l, —1,-1,-1). Furthermore, the vector e A can be ana- 
lyzed with the use of its components e^ in a coordinate basis, that 
is e A = e^. 

In such a construction, the metric tensor is obtained from the 
dual vierbein as 

g ll r(x) = JlAB<$(x)eS(x). (1) 

Contrary to General Relativity, which uses the torsion-less Levi- 
Civita connection, in the present formalism ones uses the curvature- 
less Weitzenbock connection T^ = e A d,e^ (Weitzenbock d 19231) ). 
whose torsion tensor reads 



T, 
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Moreover, the contorsion tensor, which equals to the difference 
between Weitzenbock and Levi-Civita connections, is defined as 
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Using these quantities one can define the teleparallel La- 
;rangian, which is the torsion scalar, as (Hayashi & Shirafuji 
1979J) . Maluf d 19941) . Arcos & Pereira J2004h ) 
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In summary, in the present formalism all the information concern- 
ing the gravitational field is included in the torsion tensor T A , and 



the torsion scalar T arises from it in a similar way as the curvature 
scalar arises from the curvature (Riemann) tensor. 

While in the teleparallel equivalent of General Relativity 
(TEGR) the action is just T, the idea of f(T) gravity is to gen- 
eralize T to a function T + f(T), which is similar in spirit to the 
generalization of the Ricci scalar R in the Einstein-Hilbert action to 
a function f(R). In particular, the action in a universe governed by 
f(T) gravity reads: 



/ = -i- f d 4 xe [T + f(T) - 2A + L m ] . 

ItCTlj J 



(4) 



where e = det(e^) = y[^g, G is the Newton's constant (we also 
set for convenience the light speed to one), L„, stands for the matter 
Lagrangian, and we have added for completeness a cosmological 
constant A (which could alternatively be absorbed into f{T)). We 
mention here that since the Ricci scalar R and the torsion scalar 
T differ only by a total derivative (Weinberg (2008)), in the case 
where f(T) is zero the action $4$ is equivalent to General Relativity 
with a cosmological constant. 

Variation of the action l(4) with respect to the vierbein gives 
the equations of motion 



-'d^S/OU + f, T ] - e A A T% A s;» + e p A S/ v d fl {T)f iTT 



--e A [T + f(T)-2A] 
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(5) 



where f T and f TT denote respectively the first and second deriva- 
tives of the function f(T) with respect to T. Note that the tensor 
jv\em\ on me right-hand side is the usual energy-momentum ten- 
sor. 



3 SPHERICALLY SYMMETRIC SOLUTIONS 

We are interesting in extracting spherically symmetric solutions of 
the four-dimensional f(T) gravity presented in the previous sec- 
tion, and in particular to extract the correction to the Schwarzschild 
solution of General Relativity. Up to now in the litera ture it has 
been done in three dimensions (Gonzalez et al. ( 120 121) ) and par- 
tially in four dimensio ns for s pecific f(T) ansat zes (D aouda et al. 
d201lh. Da ouda et al. d2012al) . Boehmer et al. d201ll) . Daouda et 
al. (2012b)), while there is also a different approach, namely to find 
the general (non-diagonal) vierbein choice that cor respond s exactly 
to the Sch warzsc hild solution (Ferraro & Fio rini 1 201 lot) , Ferraro 
& Fiorini d201ld) , Wang d201ll) . Maluf et al. d2012h ). 

Let us consider as usual the general metric ansatz of the form 



ds 2 = N(rfdf - K(r)~ 2 dr 2 



■ R(r) 2 dQ 2 , 



(6) 



where dQ. 2 = d6 2 + sin" 9 dip 2 and where N(r), K(r) and R(r) are 
three unknown functions. Note th at we d o not restrict in the case 
N(r) = K(r) as in (Daouda et al. ((2012bJ)), since we desire to be 
as general as possible. As a next step we need to extract the cor- 
responding vierbeins that give rise to the above metric, through re- 
lation l[T). The simplest choice is to use a diagonal vierbein of the 
form 



diag (N(r), K(r)-\ R(r), R(r) sin e) . 



(7) 



However, we mention that although in the case of linear-in-r grav- 
ity, such a simple relation between the metric and the vierbeins 
is always allowed, in the general f(T) gravity this is not the case 
anymore, and in general one has a more complicated relation con- 
necting the vierbein tetrad with the metric, with the fo rmer b eing 
non-diagonal even for a diagonal metric (Sotiriou et al. d201ll) ). as 
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it arises from Lorentz transformations of ((7). Nevertheless, in the 
4D cosmological investigations of f(T) gravity (Ferraro & Fior - 
ini J2007h. Fe rraro & Fiorini (2008), Bengochea & F erraro ( J2009h . 
Lin der J20irJ) , Mvrzakul ov J201 lh . Yerzhan ov et al J2010l) , Wu & 
Yu (1201 Pal) . Chen e t al. J201 ll) . Wu & Yu fcOllah . Bamba et al. 
201Ch . Dent et al. J2011l) ,"zheng & Hu ang 



2011 



Bamba et al. 



Yang l l201ll) . Z hang e t al. J201lh . Cai e tal. fcOllh . Chat- 
topadhy ay & Debnath feoilh. Sh arif & Ra ni (120111) . Wei et al. 
( 1201 lap. Fer raro & Fiorini J201 la), Wei et al. fcOllbl), Capo zziello 
et al. J201ll) , Wu & Yu feOllbh . Ba mba & Geng J201lh, Gen g 
et al. J201 ll) . Wei fe012h . Geng et al J2012h . Wu & Geng fcoilh . 
Bohme r et al J2012h . Karami & Abdo lmaleki fcOllh. At azadeh & 
Darabi J2012h . Farajollahi et al. J2012h . Yang et al J20121) . Karami 



& Abd olmaleki J2012ah. Kar ami & Abdolmaleki ( 




2012b), Xu et al 
20121). Liu etal. 



20121) . Bamba et a l. J2012rJ) . Setare & Hound jo _ 

Wu & Yu J2010bl). B engochea J201ll). Fer raro & Fiorini 
Daouda et al. ( 1201 ll) . Daoud a et al. J2012a|)). as well as in 



its bla ck hole solutions ( Wang J201 lh . Miao et al. d201lh . Wei et al 
I2OI2J), Ferraro & Fiorini J2011br) ). the authors still use the simple 
relation between the vierbeins and the metric, as a first approach 
to reveal the structure and the features of the theory. Therefore, in 
the present work, we also assume the simple relation between the 
vierbeins and the metric as a first approach on this novel theory, 
capable of revealing the main features of the solutions. Clearly a 
detailed investigation of the general vierbein choice in f(T) grav- 
ity, and its relation to extra degrees of freedom, is a necessary step 
for the understanding of this new theory (Cai et al. (2012)). 

Inserting the vierbein choice into the torsion scalar defini- 
tion l[3} we obtain 
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where primes denote derivatives with respect to r. Thus, substi- 
tution in the field equations {3J, neglecting the matter energy- 
momentum tensor, provides the following separate equations of 
motion: 
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Finally, we mention here that while in teleparallel gravity, as well 
as in General Relativity, the off-diagonal components of the field 
equations vanish, in f(T) gravity even the diagonal vierbein field 
gives rise to an extra equation (Boehmer et al. ( 120111) ). In particular, 
the r-6 equation reads: 



Kcote, 
2R 2 



-T'f r 



0, 



(12) 



The above equations J8b-112t cannot obtain analytical solu- 
tions in general. Thus, we proceed making two assumptions. The 
fist is the usual one in order to reduce the unknown functions, 
namely without loss of generality we set R(r) = r (going to more 



general R(r) will only lead to more complicated expressions but 
qualitatively the results will be similar). The second assumption is 
related to the f(T) ansatz. In particular, since the usual General 
Relativity is obtained by setting f(T) to zero, it is expected that a 
realistic f(T) must be small comparing to T, and as we mentioned 
this has been verified using cosmological observations (Wu & Yu 
( 1201 Obi) . Bengochea ( 120110 ). Therefore, by expanding every f(T) 
ansatz in T-powers, we deduce that as a first approximation we can 
consider a form of 



f(T) = aT 2 + 0(T 3 ), 



(13) 



where we do not write the constant and the linear terms in the 
above expansion since they can be respectively absorbed in A and 
T terms of the action (0. In the above expression a is a pure f(T)- 
parameter that determines at first approximation the divergence 
from teleparallel gravity, that is from General Relativity. We men- 
tion her e that in Friedmann-Robertson- Walker metric T = -6H 2 
(Linder yOlOh), while in the Schwarzschild metric T oc A fo r large 
r (Aldrovandi & Pereira ( |2010|) . Ferraro & Fiorini J2011bh ). that 
is in both cases T <K 1 in S.I units, which is an additional indi- 
cation that the above expansion is justified (moreover that is why 
we do not consider ansatzes of negative power laws of T). Finally, 
the above ansatz has been used as a first non-linear correc tion o f 
f(T) gravity i n (Ferra ro & Fiorini ( 1201 lbl) . Gonzalez et al. ( 120121) . 
Daouda et al. ( I2012bl) ). In summary, going beyond the assumption 
( I13t . although it would bring significant mathematical complica- 
tions, it would not qualitatively change our results. 

Thus, under the ansatz dl3t . equations ll9ll-lll2t up to O (a/r 2 ) 
lead to the solution: 



N(r)~ 



K(r) 1 
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(15) 



where we have restored the light speed, and where we have intro- 
duced the integration constant GM, with M the mass of the spher- 
ical object. We mention here that in the above solutions we have 
kept only up to linear terms in a/r 2 , since higher-order terms are 
expected to be extremely small for r in the region of the solar sys- 
tem that we have orbital data. In any case, this assumption will 
be checked a posteriori, after we extract the bounds on a. Finally, 
we straightforwardly observe that in the limit a — > 0, the above 
solution coincides with the Schwarzschild-de Sitter one (Rindler 
(2006)) as expected, while for a, A — > it becomes the usual 
Schwarzschild one. 



4 SOLAR SYSTEM CONSTRAINTS 

In the previous section we extracted the spherically symmetric so- 
lutions of f(T) gravity, in the case where there is a small but general 
divergence from the T-gravity, that is from the teleparallel equiva- 
lent of General Relativity (TEGR). Therefore, in the present section 
we apply these solutions to the gravitational field of the Sun, and 
we use t he late st solar system data from planetary orbital motions 
(Fienga ( 120 111) ) to perform a sensitivity analysis on the two model 
parameters, namely the cosmological constant A and the parameter 
a that determines the divergence from General Relativity. Finally, 
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as usual we neglect the effects of spatial curvature, which is a robust 
approximation for solar system scales. 

As usual, for a gravitational field determined by a spherically 
symmetric metric g^,, the Newtonian potential and its corrections 
are extracted from the 00 component, namely if goo = 1 +^oo> where 
Iiqq is dimensionless, then the potential is U = h^c 1 j2. Therefore, 
from the g o form of f(T) gravity spherical solutions, that is from 
N(r) 2 of d!4t . we observe that there are two corrections to the New- 
tonian potential £/n affecting a test particle's orbit, namely 



2 ,.2 



u A 



u a 



Ac ■ r 
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iac 

r 2 ' 



(16) 



(17) 



We mention here that the term -6a A is uniform, that is it does not 
depend on r, and therefore it does not affect orbital motions, while 
the term oc aAGM/r does not yield secular orbital precessions, 
since it just corresponds to a Newtonian term with GM rescaled. 
Thus, in the following it is adequate to consider only the terms 
d!6t,dl7ll. Lastly, note that the dimensions of the two parameters 



[A] 



L 2 . 



(18) 
(19) 



The orbital effects of extra-potentials having the functional 
forms of d!6t. dl7t have been incorporated analytically several 
times with a variety of different approaches in the framewor k of so - 
lar system investigations (Islam (1983), Cardona & Tejeiro ( 1998), 
Kerr e t al. |2003), Kraniotis & Whitehouse J2003h , Jetz er & S ereno 
(2006), Ka graman ova et al. d2006l). Ser eno & Jetzer J2006I), Ad - 



kins e t al. (|2007bl) , Adkin s et al. d2007<J) . Sereno & Jetzer b007l) . 
Iorio d2006h . Iorio J2008I) . Iorio J2012I) ) . and they can be straight- 
forwardly handled using, for inst ance, t he standard Lagrange per- 
turbative scheme (Bertotti et al. (2003)). From the equations for 
the variations of the osculating Keplerian orbital elements (Bertotti 
et al. (2003)) we immediately deduce that only the longitude of 
the pericenter rn = Q. + a> (with CI the longitude of the ascend- 
ing node and a> the argument of pericenter) and the mean anomaly 
M undergo secular precessions, due to the spherical symmetry of 
d!6t.dl7b and their time-independence. On the other hand, the semi- 
major axis a, the eccentricity e, the inclination / of the orbital plane 
to the reference \x,y) plane chosen, and the longitude of the as- 
cending node il, remain unaffected [J. However, from the point of 
view of a comparison with solar system observations, we are inter- 
ested only in the secular precessions of the perihelia, since the rates 
of the mean anomaly are dominated by the further, additive mis- 
modeling in the Keplerian mean motions n b — yGMa^ due to the 
uncertainty cr GM = 1 km 3 s~ 2 in the solar gravitational parameter 
GM (Konopliv et al. J201lh ). 

The Lagrange equation for the secula r precession of the lon- 
gitude of pericenter is (Bertotti et al. ( 120030 '): 



I d-nr 
\~dT 



n^w- 
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1/2 



8(H) tan(//2) d<7?> 
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(20) 



1 Indeed, the Lagrange rate equation for a (Bertotti et al. 420031) ) contains 
the partial derivative of the averaged perturbing potential (Upert) with re- 
spect to M, which is proportional to t. The Lagrange equations for e, /, H 
(Bertotti et al. (2003J)) are formed with the partial derivatives of ( J/ pcl t ) with 
respect to fl, ui, I (and Al as well in the Lagrange equation for e), which are 
absent in spherically symmetric perturbing potentials. 



Table 1. Supplementary rates of change Afl and Act, of nodes fi and perihe- 
lia w respectively, of some planets of the solar system, in milliarcseconds 
per century (mas cty~'), estimated by Fienga et al. i201lb with the IN- 
POPlOa ephemerides. Data from Messenger and Cassini were included for 
Mercury and Saturn. The reference {x, y) plane is the mean Earth's equator 
at J2000.0. 





A "(w) 


Acr(=) 


Mercury 


1.4 ± 1.8 


0.4 ±0.6 


Venus 


0.2 ± 1.5 


0.2 ± 1.5 


Earth 


0.0 ± 0.9 


-0.2 ±0.9 


Mars 


-0.05 ±0.13 


-0.04 ±0.15 


Jupiter 


-40 ± 42 


-41 ± 42 


Saturn 


-0.1 ±0.4 


0.15 ±0.65 



where 

K = (U pm ) 



(21) 



is the average of the perturbing potential over one full orbital revo- 
lution of the test particle. Inserting dl6U 17t into ( 120b we obtain 



CT A 



2\ Mb ) 



3ac 



a 4 n b (l-e 2 )' 



(22) 



(23) 



where from now on for conciseness we neglect the brackets (...) 
denoting the average over one orbital period, while the overdots 
indicate the time derivative. 

In order to preliminaril}! 2 ] infer upper bounds on both A and 
a, we compare d22t.d23t for some planets to the supplementary 
rates of change Act of their perihelia et, listed in Table Q] deter- 
mined from observations by Fienga et al. d201ll) with the INPOlOa 
ephemerides. 

However, one must be careful in straightforwardly compar- 
ing d22M23l) to the extra-rates of Table [T] since Act, by construc- 
tion, account for all those competing forces which were not mod- 
eleqjin the software of the INPOPlOa ephemerides. Moreover, Act 
are also affected by the unavoidable mismodeling in the effects 
which were actually included in the models fitted to the observa- 
tions in INPOPlOa. Thus, one should correct Act for those unmod- 
eled/mismodeled dynamical features of motion, whose magnitude 
is at least comparable to the uncertainties in determining Act them- 
selves. In this respect, the actions of the Sun's quadrupole mass 
moment J? (Rozelot & Damiani d201lh ), currently known with a 
10% uncertainty, and the general relativistic Lense-Thirring effect 
(Lense & Thirring dl918l) ), caused by the Sun's angular momentum 
5 and not modeled at all in INPOPlOa, are relevant since they in- 
duce competing secular perihelion precessions whose size is large 

2 An alternative approach would consist of explicitly modeling f(T) grav- 
ity in the planetary data processing softwares of the ephemerides, and fitting 
such ad-hoc modified dynamical models to the same observations. Thus, a 
and A would be estimated as solve-for parameters. 

3 We only deal with well known and established forces, both Newtonian 
and Einsteinian, since it would be pointless to invoke the potentially dis- 
turbing action of any sort of putative, exotic effect, whose alleged biasing 
action may, after all, be present even if a full covariance analysis is per- 
formed, as outlined in footnote [2] 
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enough, as far as the inner planets, whose extra-rates are more ac- 
curately determined, are concerned. The explicit analytical expres- 
sions, valid for a generic spatial orientation of the Sun's rotational 
axis (actually it does not point exactly to the North Celestial Pole), 
of the perihelion precessio ns due to J% and the Lense-Thirring ef- 
fect, can be found in Iorio J2011I) . 

Having at our disposal the extra-precessions of the perihelia 
of more than one planet, it is possible to set up linear combinations 
of the form 



A 



,>) 



/.<•': 



Mercury, Venus, . . ., 
;,(■■■) 



<LT) 



(24) 



where the coefficients kr come from the analytical expressions 
of the precessions due to the various effects considered simi- 
lar to J22M23t. In summary, it is possible to construct a non- 
homogeneous linear system of four equations for the four un- 
knowns A,a,J 2 ,S, which is solvable. The resulting values of A 
and a are, by construction, independent of both Ji and 5 , what- 
ever their values may be. By propagating the uncertainties in Attr A 
entering the expressions of A and a, it is possible to preliminarily 
constrain them. By using the first four inner planets, it turns out that 



|A] 



< 1.8 x 10 4 m 2 . 



(25) 



(26) 



It may b e notic ed that the bound on A of d25t is tighter than the one 
in Iorio ( 120061) by about one order of magnitude. 

Since we have now extracted the upper bound on the f(T)- 
correction allowed from solar system orbital motions, we can ex- 
amine the accuracy of the approximations made in our analysis. 
Firstly, taking r to be the Mercury mean heliocentric distance 

^Merc = «(l + 2) = ^.9 x 10 10 m ( a i s i ts semimajor axis and e 
its eccentricity), which is the smallest distance where we have ob- 
servational data, we obtain a/r^^ =5x 10~ 18 , and this justifies 

our approximation to neglect orders of 0(a/r 2 \ in d 14fc.d 15b. Ad- 
ditionally, given the small value of the cosmological constant, we 
indeed see that the two cross-terms of <14t . which are proportional 
to oA, that we did not take into account in the analysis, are many 
orders of magnitude smaller that the terms J16b.J17t that we did 
use (whether we use Mercury's data, that is enhancing the a-term, 
or the Pluto's data, that is enhancing the A-term), and thus our ap- 
proximation is well-justified. Moreover, calculating T from <[8j we 
find that T as -2A(1 + 6aA) for r » 1 (in S.I units), a result that 
verifies our indication that T ■« 1 (in S.I units) in d 1 3b - Finally, we 
can also verify that aT 2 <k T, that is the correction to the linear 
term is very small as expected. 

Relation d26t is the main result of the present work, that is 
it is the upper bound of the /(r)-correction to teleparallel grav- 
ity (that is to General Relativity), allowed from solar system or- 
bital motions. In order to provide this allowed correction in a more 
transparent way, we denote by At//(r) the divergence of the gravi- 
tational potential in f(T) gravity from that in the teleparallel equiv- 
alent of General Relativity, defined as the ratio of the correction 
term 6|a|r~ 2 in l !14t to the standard term 2GMc~ 2 r~' . The maxi- 
mal allowed value of At/^rj will be acquired at the Mercury's orbit 
T = ^Merc, in which the correction term is the largest one, since ?Merc 
is the smallest distance in the region where data exist. Thus, we find 



At/ 



/m 



6\a\r N 



2GMc- 2 T' 



<6.2xl0" 10 ; 



(27) 



as expected the allowed /(r)-correction affecting solar system or- 
bital motions is very small. 



5 CONCLUSIONS 

In this work we used latest data from solar system orbital motions 
to constrain f(T) gravity. In particular, considering the basic and 
usual ansatz f(T) = aT 2 , which is a good approximation in all re- 
alistic cases, and including also a cosmological constant A for com- 
pleteness, we extracted the spherically symmetric solutions of the 
theory, which coincide with the Schwarzschild-de Sitter one in the 
limit a — > 0. Thus, by describing the Sun's exterior gravitational 
field by these solutions, we were able to use data from planetary 
motions in order to constrain a and A. 

Concerning the cosmological constant A, we obtained the 
usual tiny bounds. Interestingly enough, our current /(r)-analysis 
leads to o ne ord er of magnitude tighter constraints than General 
Relativity (2006), that is without considering the o--term in the met- 
ric. Concerning the pure /(T)-parameter a, the obtained bound j26t 
leads to a maximal allowed divergence of the gravitational potential 
in f(T) gravity from that in the teleparallel equivalent of General 
Relativity of the order of < 6.2 x 10~ 10 for the smallest distance 
'Wc (Mercury mean heliocentric distance) that we have data, and 
where the correction term is the largest one. 

The above obtained small divergence from General Relativity 
is much smaller than the corresponding (significantly small too) di- 
vergence that is predicted from cosmological constrains (Wu & Yu 
(2010b), Bengochea (2011)), where one also finds that f(T) must 
be close to the linear-in-r form. This is standard in observational 
constraints and justified, since the solar system observations are al- 
ways more accurate than the cosmological ones. In summary, in the 
present analysis we did verify, in a different context, the expected 
result that the allowed divergences of f(T) gravity from the linear 
(teleparallel) form are significantly small, and this should be taken 
into account in f(T) model building. 

In the above analysis we remained in the diagonal vierbein 
ansatz, since in this case one can safely elaborate the spherical so- 
lutions. However, observing the features of some specific spherical 
solutions under non-diagonal ansatzes (Daouda et al. ( I2012bh ). in 
which one obtains similar terms with the present expressions, we 
deduce that the above analysis would lead to qualitatively similar 
results even in those cases. However, since for the moment it is 
not clear how to elaborate the non-diagonal cases and in particular 
how to handle the extra degrees of freedom at the background and 
especially at the perturbation level, we preferred to remain in the 
diagonal vierbein scenario. 

We close this work by mentioning that since in the solutions 
of the aforementioned analysis there appear terms of the form a/r 2 
(which are small for planetary motions as discussed above a nd eve n 
a full Parametrized-Post-Newtonian (PPN) analysis (Will J2005h ) 
will not change the obtained results), they could be significant 
at much smaller distances, outside the applicability region of the 
present analysis. Thus, an interesting and necessary investigation 
would be to constrain the allowed correction of /(r)-gravity using 
different scenarios like, for example, fast extrasolar planets orbiting 
their parent stars at distances smaller than Mercury. Furthermore, 
Earth-based laboratory experiments (in which the simple spherical 
geometry will not be the only case), where the allowed corrections 
are also expected to be small according to the recen t tests on Gen- 
eral Relativity (Dimopoulos et al. J2007I) . Turyshev J2008I) ). would 
be worthwhile. However, what could happen in even smaller dis- 
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tances is unknown, where the possibility of large divergences from 
General Relativity could remain open. 
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